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Abstract We report a detailed investigation of the bulk motions of the nearby Galactic 
stellar disk, based on three samples selected from the LSS-GAC DR2: a global sample 
containing 0.57 million FGK dwarfs out to ~ 2 kpc, a local subset of the global sample 
consisting ~ 5,400 stars within 150 pc, and an anti-center sample containing ~ 4,400 
AFGK dwarfs and red clump stars within windows of a few degree wide centered on the 
Galactic anti-center. The global sample is used to construct a three-dimensional map of 
bulk motions of the Galactic disk from the solar vicinity out to ~ 2 kpc with a spatial 
resolution of ~ 250 pc. Typical values of the radial and vertical components of bulk 
motion range from —15 km s” 1 to 15 km s _1 , while the lag behind the circular speed 
dominates the azimuthal component by up to ~ 15 km s _1 . The map reveals spatially 
coherent, kpc-scale stellar flows in the disk, with typical velocities of a few tens km s _1 . 
Bending- and breathing-mode perturbations are clearly visible, and vary smoothly across 
the disk plane. Our data also reveal higher-order perturbations, such as breaks and ripples, 
in the profiles of vertical motion versus height. From the local sample, we find that stars 
of different populations exhibit very different patterns of bulk motion. Finally, the anti- 
center sample reveals a number of peaks in stellar number density in the line-of-sight 
velocity versus distance distribution, with the nearer ones apparently related to the known 
moving groups. The “velocity bifurcation” reported by Liu et al. (2012) at Galactocentric 
radii 10-11 kpc is confirmed. However, just beyond this distance, our data also reveal a 
new triple-peaked structure. 
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1 INTRODUCTION 

The stellar velocity distributions encode the assemblage history and structural information of the Milky 
Way. The Milky Way is often approximately considered to be steady and axisymmetric. In such a model, 
disk stars are expected to have quite smooth and single-peaked velocity distributions, without significant 
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bulk motions in either the radial or vertical directions (see e.g. Schwarzschild & Villiger 1907; Schonrich 
& Binney 2012). 

This simple picture is, however, challenged by the recent detections of significant bulk motions 
in the nearby stellar disk. Using ~ 1 x 10 4 stars from the Sloan Extension for Galactic Understanding 
and Exploration survey (SEGUE; Yanny et al. 2009), Widrow et al. (2012) determine the vertical bulk 
motions as a function of height from the disk mid-plane. Their results reveal wave-like perturbations 
and some differences in bulk motion between blue and red stars. Williams et al. (2013) report the bulk 
motions of ~ 7 x 10 4 red clump stars selected from the Radial Velocity Experiment (RAVE; Steinmetz et 
al. 2006). They find a north-south asymmetry both in the radial and vertical bulk motions; and curiously, 
the vertical bulk motions exhibit some compression and rarefaction patterns across the disk. North-south 
asymmetry and compression patterns are also observed by Carlin et al. (2013) with ~ 4xl0 5 F-type 
stars selected from the LAMOST Galactic spectroscopic surveys (Cui et al. 2012; Zhao et al. 2012) 
and RAVE. Bovy et al. (2015) report fluctuations in differences between the observed and predicted 
stellar velocity distributions on scales ranging from 25 pc to 10 kpc, using data from the Apache Point 
Observatory Galactic Evolution Experiment (APOGEE; Eisenstein et al. 2011), RAVE and the Geneva- 
Copenhagen Survey (GCS; Nordstrom et al. 2004). They attribute the fluctuations to the effects of the 
central bar. Besides, Liu et al. (2012) reveal a feature of “velocity bifurcation” at Galactocentric radii 
10-11 kpc with ~ 700 red clump stars in the Galactic anti-center direction, and interpret it as a resonance 
feature of the central bar. Theoretically, Siebert et al. (2011), Debattista (2014) and Faure et al. (2014) 
show that bulk motions can be excited by spiral arms, and Widrow et al. (2014) investigate the dynamical 
effects of a passing satellite or dark matter subhalo. 

By June 2014, the LAMOST Spectroscopic Survey of Galactic Anti-Center (LSS-GAC; Liu et 
al. 2014; Yuan et al. 2015) has obtained spectra with S/N(4650A) >10 for over one million stars. The 
second data release of value-added catalogs (DR2) provides their line-of-sight radial velocities, stellar 
atmospheric parameters, three-dimensional extinctions, distances and proper motions. In this paper, we 
report an investigation of bulk motions based on three samples selected from the LSS-GAC DR2: A 
global, a local and an anti-center sample. The global sample, containing ~ 0.57 million FGK dwarfs, is 
used to obtain a three-dimensional map of bulk motions of the Galactic disk from the solar neighborhood 
out to ~ 2 kpc with a resolution of ~ 250 pc. The second sample, a subset of the global one, contains 
~ 5,400 stars within 150 pc; we use it to investigate the bulk motions of different stellar populations 
in the solar neighborhood. Finally, we investigate whether the “velocity bifurcation” reported by Liu et 
al. (2012) is a localized structure via an analysis of the line-of-sight radial velocities of the anti-center 
sample, consisting ~ 4,400 AFGK dwarfs and red clump giants, tracing a contiguous distance range of 
4 kpc in the Galactic anti-center direction. The paper is organized as follows. After a description of the 
data in Section [2] we present results based on the three samples in Sections [3}{5] in turn. We then close 
with a brief summary in Section[6] 

2 DATA 

2.1 LSS-GAC 

Data used in this work are from the LSS-GAC DR2 (Xiang et al. 2015b, in preparation). LSS-GAC - the 
LAMOST Spectroscopic Survey of the Galactic Anti-Center (Liu et al. 2014, Yuan et al. 2015) surveys a 
3400 deg 2 contiguous area centered on the Galactic anti-center (150° < l < 210°, —30° < b < 30°) and 
aims to obtain optical (AA 3700 — 9000), low-resolution (R ~ 1,800) spectra for a statistically complete 
sampltQ containing ~ 3 million stars of all colors and of magnitudes ranging from r = 14.0 to 17.8 mag 
(down to 18.5 mag for a limited number of fields). Over 1.5 million stars of —10° < Dec < 60° and 
brighter than r = 14 mag, are also selected with a similar selection algorithm and observed under bright 
lunar conditions. LSS-GAC also contains a sub-program, targeting various types of objects in the M31- 

1 In the sense that the sample stars are selected uniformly on celestial sphere to given limiting magnitudes and randomly in 
color-magnitude diagrams, such that with due considerations of the various well-defined selection effects, the spectroscopically 
targeted sample stars can be used to recover the underlying photometric populations (cf. Yuan et al. 2015 for detail). 
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M33 region. By June, 2014, LSS-GAC has obtained over two million spectra with a signal-to-noise ratio 
at 4650 A, S/N(4650A) > 10. 

Line-of-sight radial velocities, V r , and stellar atmospheric parameters (effective temperature T e s, 
surface gravity log g and metallicity [Fe/H]) are estimated with the LAMOST Stellar Parameter Pipeline 
at Peking University (LSP3; Xiang et al. 2015a). For S/N(4650A) > 10, LSP3 has achieved accuracies 
better than 5.0 km s -1 , 150 K, 0.25 dex and 0.15 dex for V r , T e g, log g and [Fe/H], respectively. Values 
of the interstellar extinction, with typical errors of ~ 0.04 mag, are obtained for individual stars with 
the standard pairing technique (Stecher 1965, Massa et al. 1983, Yuan et al. 2013). Distances have been 
derived based on an empirical relation between absolute magnitudes (go, ro, Jo, Ho, K s o) and stellar 
atmospheric parameters (T e g, log g, [Fe/H]). Typical distance uncertainties are 10-15% for dwarfs and 
20-30% for giants. 

Proper motions, p a cosS and //,,, are from the UCAC4 (Zacharias et al. 2013) and PPMXL (Roeser 
et al. 2010) catalogs. UCAC4 is an all-sky, astrometric catalog containing over 113 million objects 
complete to R ~ 16 mag. More than 105 million stars in UCAC4 have proper motions, with typical 
random errors of 4 mas yr _1 . Compiled based on USNO-B1.0,2MASS and a previous version of PPMX, 
PPMXL contains over 900 million sources down to V = 20 mag. Almost all of them have proper motions, 
with typical random errors ranging from 4 to more than 10 mas yr _1 . The position-dependent systematic 
errors of proper motions of UCAC4 and PPMXL have been corrected for by Huang et al. (2015) and 
Carlin et al. (2013), respectively, using ~ 1,700 or ~ 10 5 quasars as anchors. After the corrections, 
however, some unaccounted for color- and magnitude-dependent systematic errors remain, which are 
further discussed in Section [3)21 

2.2 Coordinate Systems and Galactic Parameters 

We use two sets of coordinate systems in this paper: (a) A right-handed Cartesian coordinate system ( X, 
Y, Z) centered on the Sun, with X increasing towards the Galactic center, Y in the direction of Galactic 
rotation and Z the height from the disk mid-plane, positive towards the north Galactic pole; and (b) A 
Galactocentric cylindrical coordinate system (R, <t>, Z), with R the Galactocentric distance, $ increasing 
in the direction of Galactic rotation and Z the same as that in the Cartesian system. We assume that the 
Sun has a Calactocentric distance Rq = 8 kpc. This places the Galactic center at (Xqc, Ycc, -Zgc) = (8, 
0, 0) kpc and the Sun at (Rq, 4>q, Zq) = (8 kpc, 0°, 0 kpc). The three velocity components are denoted 
as (U, V, W) in the Cartesian system and (Vr, Ivj, , Vz) in the cylindrical system. We use a flat rotation 
curve of 14 = 220 km s -1 , and solar motions (Uq, Vq, Wq) = (7.01, 10.13, 4.95) km s -1 (Huang et 
al. 2015) relative to the Local Standard of Rest (LSR). 

2.3 Samples 

2.3.1 The global sample 

We use a global sample of FGK dwarfs to construct a three-dimensional map of bulk motions of the 
Galactic disk from the solar vicinity out to ~ 2 kpc. The stars are selected from the LSS-GAC DR2 
by applying effective temperature and surface gravity cuts, 4,200 < T e g < 6,800 K and 3.8 < log g 
< 5.0 dex. Stars without UCAC4 or PPMXL proper motions are excluded and duplicate ones of lower 
spectral SNRs removed. This leads to a sample of ~ 0.70 million FGK dwarfs. 

We derive the three-dimensional positions, (X, Y, Z) and (R, 't>, Z), and three-dimensional veloc¬ 
ities, (U, V, W) and (Vr, V$, Vz), of individual stars from their celestial coordinates, (l, b), distances, 
d, line-of-sight radial velocities, V r , and proper motions, g a cosS and gs- For proper motions, both 
UCAC4 and PPMXL catalogs are used, yielding two independent sets of velocities for each star. We use 
a Monte Carlo method of error propagation to estimate the random errors of positions and velocities. 
This is done by randomly sampling the distance moduli, radial velocities and proper motions, assum¬ 
ing Gaussian random error distributions. For each realization of a given star, we calculate the position 
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Fig-1 Color-magnitude diagram of the global sample (contours), and the median differences 
of UCAC4 and PPMXL proper motions (//,, cos(<5) on the left and ps on the right) of the 
member stars falling in the individual color-magnitude bins (colorscale). The bins have a 
width of 0.1 mag in g — r and 0.5 mag in r. The contours increase from 100 to ~ 1 x 10 4 stars 
per bin with equal logarithmic steps. The colorbar is in units of mas yr -1 . 


and velocity. This is repeated 1000 times to obtain their distributions. The standard deviations of the 
distributions are adopted as random errors of the position and velocity of the star in concern. 

We then remove outliers and stars with large random errors by requiring: —200 <Vr < 200 km s _1 , 
0 < V,p < 400 km s _1 , —200 < Vz < 200 km s _1 and random errors in Vr, V$ and Vz no larger than 
50 km s _1 . Velocities obtained with the UCAC4 and PPMXL proper motions should both meet these 
criteria to retain a star in the sample. This leaves ~ 0.57 million FGK dwarfs in the final sample. The 
contours in Fig. Q] display the sample in the (g — r, r) color-magnitude diagram. The sample covers a 
range of 10-17 mag in r-band magnitude and 0.1-1.7 mag in g — r color. The spatial distributions of 
the stars can be found in Fig. [2] which shows that the stars occupy a significant volume of the Galactic 
disk, spanning from the solar vicinity out to ~ 2 kpc. The global sample is three-dimensionally sliced 
into small bins, with binsize (A R, A<f>, A Z) = (250 pc, 2°, 250 pc). The entire volume of Fig. [2] thus 
contains a total of 2304 bins. Amongst them, 662 bins are found to contain stars, ranging from several 
to a few 10 4 ; and 252 bins have no fewer than 100 stars. 

2.3.2 The local sample 

A local sample of stars are selected by applying a distance cut, d < 150 pc, to the global sample. The 
local sample is less affected by the potential systematic errors in proper motions due to the distance 
cut. The sample contains 5,443 FGK dwarfs, and is further divided into six populations according to 
their effective temperatures (F/G/K types) and metallicities (metal-rich and metal-poor). Table 1 (see 
Section® lists the numbers of stars, along with the effective temperature and metallicity ranges of each 
population. 

2.3.3 The anti-center sample 

The anti-center sample is consisted of four populations, the GK-, F-, A-type dwarfs and red clump (RC) 
giants. The dwarfs are selected from the LSS-GAC DR2 with the effective temperature and surface 
gravity cuts, 4,200 < T e g < 6,000 K and 3.8 dex < log g < 5.0 dex for GK-type stars, 6,000 < T e g 






Fig. 2 Spatial distributions of the global sample. Greyscale shows the logarithmic surface 
density of stars in units of pc” 2 (left and middle panels) or in pc" 1 deg” 1 (right panel). With 
spacings (A R, A<f>, A Z) = (250 pc, 2°, 250 pc), the magenta grids show how the sample is 
divided three-dimensionally when we construct the bulk motion maps. 


< 6,800 K and 3.8 < log g < 5.0 dex for F-type stars, and 6,800 < T e g < 10,000 K and 3.2 < log g 

< 5.0 dex for A-type stars. The RC giants are selected by Huang et al. (this volume), also from the 
LSS-GAC DR2, based on their positions in the [Fe/H]-T e ff-log g and [Fe/H]-color parameter spaces. 
The PARSEC stellar evolution models (Bressan et al. 2012) and high-quality asteroseismology data of 
Kepler are used to define the location of RC giants in the parameter spaces. The contamination and 
completeness of the selection algorithm are less than 5% and greater than 90%, respectively. RC stars 
thus selected have distance accuracies better than 5-10%. From the chosen AFGK-type dwarfs and 
RC giants, we then select those falling in windows centered on the Galactic anti-center, \l— 180°|< <5, 
\h\< 5, where <5 is set to 4° for GK- and F-type dwarfs, 6° for A-type ones and 2° for RC giants, a 
compromise between a sufficient number of stars (preferring a large 5) and a good alignment between 
the Galactocentric and line-of-sight radial velocities (preferring a small 5). The final anti-center sample 
contains 1,716 GK-, 1,402 F-, 359 A-type dwarfs and 958 RC giants. With increasing luminosities, these 
four populations of stars are able to trace the Galactic disk in the anti-center direction from the solar 
vicinity out to a Galactocentric distance of 12 kpc. The distance distributions of the anti-center sample 
can be found in Fig.fl7l(Sectionl5Ti. 

3 THE THREE-DIMENSIONAL BULK MOTION MAP 
3.1 Characterizing Bulk Motions 

To construct the bulk motion maps, we three-dimensionally slice the global sample into small bins (see 
Section l2.3.1l and Fig. O. Fig.[3]displays the velocity distributions of two typical bins. The first bin is a 
close one and thus less affected by the random errors propagated from the proper motions and distances. 
As Panel (a) shows, the distribution of Vr in the first bin shows deviations from a simple, smooth 
Gaussian distribution, including two peaks at —15 km s” 1 and 10 km s” 1 , respectively. The peaks are 
produced by the known moving groups (see e.g. Dehnen 1998; Antoja et al. 2012): the —15 km s” 1 
peak comes from the Sirius and UMa groups, while the one at 10 km s” 1 from the Pleiades and Coma 
Berenices groups. Deviations such as those are however not so big - a single Gaussian can still fit well 
the overall distribution and the median Vr finds the Gaussian peak quite precisely. A visual examination 
of the 252 bins with no fewer than 100 stars has found clear signatures of moving groups in the Vr 
profiles of only 4 bins, all near the Sun. Again, the presence of those moving groups hardly affect the 
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determinations of the centroids of the distributions through the median values of Vr. Those moving 
groups, however, do not show up in Vj, as obviously as in Vr of the same bin [Panel (b)]. They are 
probably lost in the much narrower V$ distribution (The moving groups have typical widths in Vr or 
Vj, of > 10 km s” 1 ; see e.g. Fig. 3 of Dehnen 1998). Unlike Vr, has long been known to skew 
towards lower values. This is clearly seen in the current case. Still, a single Gaussian fits the data and 
the median value finds the peak position well. Panel (c) shows quite a simple, single-peaked, symmetric 
profile of Vz of the first bin, whose peak can be easily found by Guassian fitting or by the median value 
despite the weak, extended wings at \Vz\ > 30 km s” 1 . In the second bin of farther distance [Panels 
(d-f)], large random errors propagated from the proper motions and distances would have blurred the 
signatures of any moving groups, the skewness or the extended wings in the distributions of Vr, V$ and 
Vz, respectively. In this case, the distribution profiles are well fitted by a Gaussian, whose peaks coincide 
well with the corresponding median values. Above all, as a first approximation, the bulk motions of stars 
in the individual bins can be characterized by the median values of the three components of velocity, 
denoted in this paper by (Vr), (V$) and (Vz)- To calculate those median values, we have required the 
bins to contain at least 100 stars, leaving 252 bins usable. The bulk motions thus derived have random 
errors less than ~ 5 km s , with typical values of < 2 km s’ 1 for (Vr) and < 3 km s” 1 for {V',j>) and 
(Vz). 


3.2 Systematic Errors 

3.2.1 Distances 

The random errors of stellar distances do not affect the bulk motions systematically. Slightly evolved 
stars in the global sample, with 3.8 < log g < 4.2 are found to have systematically underestimated 
distances by ~20% compared to those derived by isochrone fitting (Xiang M.-S. et al., private com¬ 
munications). However, the fraction of those evolved stars in our sample is small. Moreover, a ~20% 
underestimation in distances corresponds to an underestimation of the same percentage in the mean 
tangential velocities, which in turn contributes very minor to (Vr) and (Vz). 

3.2.2 Line-of-Sight Radial Velocities 

A comparison of the line-of-sight radial velocities of common stars observed by LAMOST and 
APOGEE has shown that the LAMOST wavelength calibration suffers from glitches of the order of 
~ 2 km s” 1 on a timescale of several days (Xiang M.-S. et al., private communications). The sample 
stars inside each bin were observed on nights ranging from typically ~ 5 up to ^ 60 days, with typical 
timespans ranging from 1-3 months to 1-3 years. We expect that the effects of those glitches largely 
cancel out on such timescale, and thus hardly affect our estimates of the bulk motions. Unresolved binary 
stars can potentially affect the line-of-sight radial velocities due to the orbital velocities of the binary 
stars around their centers of mass. However, this effect is also unlikely to affect the determinations of 
the median velocities (Vr), (I/j>) and (Vz) in a systematic way, under the reasonable assumption that 
the orbital orientations and phases of binary systems are randomly distributed. 

3.2.3 Proper Motions 

The largest systematic errors affecting our results come from the proper motions, more significantly 
at larger distances. The colorscale in Fig. |T| shows the median differences of the UCAC4 and PPMXL 
proper motions for the global sample. It can be seen that the differences vary systematically with g — r 
color and r magnitude. Typical differences for /i„cosc) range from — 1 to 0 mas yr” 1 except for very blue 
or bright stars ( g — r <0.3 mag or r < 12.5 mag), for which the differences amount to —2 mas yr -1 . For 
ps, the systematic differences typically range from — 1 to 1 mas yr” 1 ; but very red or faint stars ( g — r> 

1.3 mag or r > 16.0 mag) are found to have values up to 2 mas yr” 1 . Thus, even after corrected for 
systematics using samples of distant quasars, the proper motions from UCAC4 and PPMXL still suffer 
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Fig. 3 Distributions of Vr (left panels), V$ (middle panels) and Vz (right panels), derived 
using the UCAC4 proper motions in bins 8 < R < 8.25 kpc, 0<$<2°, 0 < Z < 0.25 kpc 
(upper panels) and 8.25 < R < 8.5 kpc, 4 < <f> < 6°, —0.5 < Z < —0.25 kpc (lower panels). 
In each panel, the histogram binsize is set to no smaller than the random errors of the 90th 
percentile of stars ordered by increasing error. The binsize is set in such a way to avoid 
artifacts introduced by the velocity errors of individual stars. The red curve is a Gaussian fit 
to the distribution of which the median value is marked by the red straight line. In Panel (a), 
Peak 1 at —15 km s -1 originates from the moving groups Sirius and UMa, and Peak 2 at 
10 km s -1 comes from Pleiades and Coma Berenices. 


from the unaccounted for color- and magnitude-dependent systematic errors. In the current work, we 
have used proper motions from both catalogs and deem result robust if it is seen using both UCAC4 and 
PPMXL. 

3.2.4 Algorithm Biases 

As we shall see later, (' Vr ), (V$) and (Vz) all vary spatially in three dimensions, i.e. in (R, 4>, Z ), or 
{X, Y, Z). If the bulk motions are estimated by projecting the stars to a two-dimensional plane, the 
results are unlikely to be a good representation of the distribution in case where the bulk motions vary 
significantly in the projected dimension. For example, the mean Vr of stars within bin (li.\ - f? 2 , Z\ - 
Z 2 ) is a weighted mean of (Vr)(<1>|.Ri - R 2 , Z\ - Z 2 ) with weights n(<3?|.Ri - R 2 , Z\ - Z 2 ), where n is 
the stellar number density of the sample used. Samples of different n could yield quite different results. 
This problem is naturally avoided by three-dimensional binning, as we have done here for the global 
sample. 
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3.2.5 Selection Effects 

As we shall see in Section|4j different populations seem to have different bulk motions. The differences 
in (Vr) and (Vz) are relatively small, i.e. on the order of a few km s _1 . However, (V$) has larger 
differences, ~ 10 km s _1 , between different populations. The global sample contains stars spanning 
a wide range of effective temperatures and metallicities, and thus suffer from selection effects. The 
selection effects in (Vr) and (Vz) are less important, comparing with the systematic errors introduced 
by proper motions (from several to more than 10 km s -1 ; see Figs.[4]|9]i- In contrast, the selection effects 
in (V$) are expected to be more significant. We do not correct for the selection effects in the current 
paper. 

3.2.6 The Solar Motions 

Bovy et al. (2015) point out that bulk motions have a significant effect on the determination of the 
solar motion with respect to the LSR. With stars spanning from the solar vicinity out to ~ 5 kpc, 
they determine the azimuthal component of solar motion to be Vq ~ 24 km s _1 , much larger than 
the previously locally-determined results. The discrepancy can be explained if the local Galactic disk 
has a ~ 10 km s -1 bulk motion (Bovy et al. 2015). In this paper, we have calculated stellar velocities 
adopting a locally-determined set of solar motions, (Uq, Vq, Wq) = (7.01, 10.13,4.95) km s -1 (Huang 
et al. 2015). Here we do not aim to obtain a better set of solar motions, which may well vary depending 
on the spatial extension of stars used to determine them in view of the finding of Bovy et al. In any case, 
the bulk motions, (Vr), (V ,\,) and (Vz), reported here can be updated in the future by simply applying 
the corresponding offsets, 5Uq, SVq and SWq, between the future new and the current solar motions, 
as we approximately have S((Vr) — Uq), <5((V$) — Vq) and S((Vz) — Wq) ~ 0 km s -1 for our sample. 

3.3 The Bulk Motions 

3.3.1 Radial Component (Vr) 

Figs. Q] and [5] show the maps of radial bulk motion (Vr). For most bins, the differences between the 
UCAC4 and PPMXL results are no larger than 5 km s -1 . (Vr) is clearly nonzero, but ranges from 
~ —15 km s _1 to ~ 15 km s _1 , consistent with the previous studies (e.g. Williams et al. 2013; Carlin et 
al. 2013). The majority parts of the disk mapped have outward bulk motions. However, inward motions 
dominate a large region near the Sun, i.e. the region of 8 < R < 9 kpc, 0 < <£> < 6°, —0.75 < Z < 
0.5 kpc, or roughly — 1 < X < 0 kpc, 0 < Y < 1 kpc, —0.75 < Z < 0.5 kpc. 

3.3.2 Azimuthal Component (V$) 

Figs. [6] and [7] show the maps of azimuthal bulk motion ( V ,\,}, from which we have subtracted a constant 
circular speed of V c = 220 km s -1 for clarity. The UCAC4 and PPMXL proper motions give consistent 
results for most bins, with differences no larger than 5 km s . The Galactic disk is dominated by 
negative values of (V$) — V c , which decreases systematically as one moves away from the mid-plane. 
Near Z ~ 0 kpc, the stars lag behind the circular speed by ~ 8 km s -1 ; but at \Z\ ~ 1 kpc, the stars have 
( Vq ,) — 14 < —15 km s -1 . Similar trends are also unveiled by Williams et al. (2013). As mentioned in 
Section l3.2.5l (V$) could be affected by some selection effects. The later-type stars, which are expected 
to have larger asymmetric drifts, are concentrated close to the mid-plane. Thus, (V$) at larger values of 
\Z\ could be slightly over-estimated due to the selection effects. 

3.3.3 Vertical Component (Vz) 

Figs. [8] and 0 show the maps of vertical bulk motion (Vz). Similar to (Vr), (Vz) ranges approximately 
from —15 to 15 km s -1 , again in agreement with the previous studies (Widrow et al. 2012; Williams 
et al. 2013; Carlin et al. 2013). However, the vertical bulk motions are more severely affected by the 
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Fig. 4 Distributions of (Vpt) obtained with either the UCAC4 or the PPMXL proper motions 
and their differences in units of km s -1 viewed in the X—Y plane. Boxes in black show bins 
where the differences are no larger than 5 km s _1 . Numbers in the brackets show the ranges 
of Z in units of kpc. The black dot denotes the location of the Sun. 
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systematic errors in proper motions, and there are more bins where UCAC4 and PPMXL give discrepant 
values. Williams et al. (2013) and Carlin et al. (2013) have found that the Galactic disk exterior to 
the solar position is compressing, i.e. stars above the mid-plane move downwards while those below 
upwards. Our results confirm that the disk sector of 0 < $ < 6° is indeed compressing, but that of —4 

< < 0° does not. The UCAC4 results even show that the disk sector of 4 < <1> < —2° is probably 
rarefying. Stars within 9 < R < 10 kpc and —2 < •!> < 0°, above or below the disk mid-plane regardless, 
are found to be co-moving upwards, indicative of a bending-mode perturbation. More discussions on 
the vertical perturbations can be found in Section [331 

3.4 Giant Stellar Flows in the Disk 

The bulk motion maps unveil spatially coherent, kpc-scale flows in the disk, as can be seen in Figs. [TO] 
and 1 1 1 1 At most locations, the UCAC4 and PPMXL proper motions give compatible results in terms of 
both speeds and directions. Fig.[l0]displays (Vr) and (Vz) as vectors in the R-Z plane. From $ = —8° 
to $ = —2°, a giant flow, of kpc scales and typical velocities \j(V r ) 2 + (Vz ) 2 > 10 km s -1 , streams 
from the solar neighborhood out- and up-wards. The nearer parts of the flow (R < 9 kpc) become 
weaker in terms of typical velocity in the sector —2 < $ < 0° and the direction is reversed between <b 
= 0 and 2°: the flow now runs from outside into the solar neighborhood. For <1> > 2°, we see the flow 
increases its speed back to > 10 km s _1 , reaching a maximum of ~ 20 km s -1 . 

Fig. QT| which plots (Vr) and (V$) — V c in the X—Y plane, displays the giant flow from an¬ 
other perspective. Due to the large asymmetric drifts, the giant flow lags behind the circular speed by 
> 10 km s' 1 at Z < —0.5 kpc or Z > 0.75 kpc. Near the mid-plane, one sees an incomplete, clockwise 
vortex structure around —1 < X <0 kpc and — 1 < Y < 0.5 kpc. It has a diameter of ~ 1 kpc and a 
typical velocity a/ (Vr ) 2 + ((V$) — V c ) 2 <10 km s -1 . This structure is most obvious between 0 < Z 

< 0.5 kpc. 

3.5 Vertical Perturbations 

The vertical kinematics has drawn considerable interest. Here we use (Vz)(Z\R, 4>) to denote the ver¬ 
tical bulk motion (Vz) as a function of height Z at a specific location (R, <f>) in the disk plane. It 
has been shown that (Vz)(Z\R, 4>) can be decomposed into the bending and breathing modes plus the 
higher-order modes (Mathur 1990; Weinberg 1991; Widrow et al. 2014). The bending mode describes a 
situation where stars below and above the disk mid-plane move synchronously upwards or downwards, 
whereas the breathing mode refers to a scenario where stars below and above the mid-plane move in 
opposite vertical directions. The disk is said to be in compression when stars above the mid-plane move 
downwards and those below upwards, or in rarefaction vice versa. 

To better display (Vz)(Z\R, <b), we keep the binsize in the disk plane A R = 250 pc and A<!> = 2° 
unchanged but increase the sampling frequency in height Z. Bins in Z are now spaced every 50 pc and 
have a height Z x 15% (corresponding to the random errors) but no smaller than 100 pc. As before, 
(Vz) is estimated by the median Vz of stars in each bin. The curves of (Vz)(Z\R, <!>) are displayed in 
Fig.ll2lfor cases where there are at least 10 data points and the range of heights spans at least 1 kpc. 

3.5.1 Bending and Breathing Modes 

Despite the complicated behavior of (Vz)(Z\R, 4>), we find that a simple linear fit can be applied to 
roughly describe their trends (see also Widrow et al. 2014), i.e. 

(V z )(Z | .R, <E>) = A(R, $)xZ + B(R, $) (1) 

where the intersection B, in units of km s -1 , parameterizes the bending mode (positive for upward and 
negative for downward co-movement) and the slope A, in units of km s” 1 kpc -1 , parameterizes the 
breathing mode (positive for rarefaction and negative for compression). 
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The distributions of the breathing-mode parameter A is displayed in Fig. [13] The figure once again 
confirms the compression and rarefaction patterns already described in Section [3.3.3l but denoted here by 
reddish and blueish cells. Typical values of A range from ~ —15 km s -1 kpc -1 to ~ 15 km s _1 kpc -1 
(with the UCAC4 proper motions) or to ~ 0 km s -1 kpc” 1 (with the PPMXL proper motions). Out of 26 
cells in total, 18 have the UCAC4 or PPMXL results consistent with each other within 5 km s -1 kpc -1 . 
Results from both sets of proper motions reveal a smooth, gradual change of A on kpc scales across the 
disk plane. The gradient is estimated to be roughly 18 km s - 1 kpc ~ 2 . 

Fig. [14] shows the distribution of bend-mode parameter B, also found to possess smooth spatial 
variations. In the case of UCAC4 proper motions, the Galactic disk is found to have both upward or 
downward bending modes. However, with the PPMXL proper motions, we find that almost the whole 
region investigated has an upward bending mode only, except for one cell in the solar vicinity with 
B ~ — 5 km s -1 . Still, more than two thirds of the cells have the UCAC4 or PPMXL results consistent 
with each other within 5 km s -1 . A gradient of ~ 15 km s - 1 kpc -1 is estimated for the bending mode 
parameter. 

3.5.2 Breaks and Ripples 

In addition to the bending and breathing modes, the vertical bulk motion exhibits some higher-order 
perturbations, such as breaks or ripples along the (Vz)(Z\R, $) curves. A representative example of 
the breaks can be seen in Panel (1) of Fig.[l2](8.5c < R < 8.75 kpc, 0 < 4> < 2°). (Vz) first decreases 
smoothly from Z = —0.5 kpc to Z = 0.6 kpc, after which it jumps sharply by ~ 10 km s -1 within 0.2 kpc 
and then comes back again to a smooth gradual decrease with Z at Z > 0.8 kpc. The entire curve is 
broken by the two turning points at Z = 0.6 and 0.8 kpc into three segments, each well approximated by 
a linear function. Similar breaks can also be found in Panels (a, c, f). 

The second kind of higher-order perturbations are characterized by the high-frequency ripples in 
the ( Vz)(Z\R , $) curves, especially at large distances from the Sun, e.g. Panels (o-z) of Fig. [12] The 
ripples have peak-to-peak amplitudes A pp up to ~ 10 km s -1 and frequencies kz of roughly 2-5 kpc -1 . 
Similar ripples have been found by Widrow et al. (2012), with comparable peak-to-peak amplitudes 
and frequencies. One may question whether these ripples arise from the unaccounted for color- and 
magnitude-dependent systematic errors in the proper motions. However, we believe that those ripples 
are probably real, as the color and magnitude distributions of stars vary smoothly with Z. An example 
is given in Fig. [15] which shows the distributions of g — r colors and r-band magnitudes of stars of 
9.25 < R <9.5 kpc and —2 < $ < 2° [corresponding to Panel (w) of Fig. [12) ■ It can be seen that both 
g — r and r vary smoothly with Z , with no clear features corresponding to the peaks or dips of (Vz) at 
Z = 0, 0.1, 0.3, 0.5 or 0.8 kpc. A small feature in g — r at Z = —0.2 kpc is possibly related to a small 
dip at the same location in (Vz) obtained with the UCAC4 proper motions. Besides, the locus of g — r 
and r median values in the (g — r, r ) Hess diagram occupies a region where the differences of p a cos6 
between UCAC4 and PPMXL are almost constant while the differences of ps vary smoothly from 3 to 
—2 mas yr -1 . Although we do not know the systematic errors of UCAC4 or PPMXL proper motions 
themselves (we only know their differences), it is unlikely that they all vary quickly along the locus while 
at the same time generating smoothly varying differences. Thus it seems to us that the systematic errors 
of proper motions may not be responsible for the observed high-frequency ripples along the (Vz)(Z\R, 
<I>) curves. 

4 BULK MOTIONS OF THE LOCAL STELLAR POPULATIONS 

It is an interesting problem whether different stellar populations share the same bulk motions. Smith et 
al. (2012) and Widrow et al. (2012) show that there are differences of up to ~ 5 km s -1 between stars of 
different metallicities or colors. Such a study is not feasible for the global sample, as the potential color- 
dependent systematic errors in proper motions may generate false bulk motion differences amongst 
populations of different colors, especially at large distances. We thus restrict ourselves to the local 
sample, which contains 5,443 FGK stars of distances < 150 pc from the Sun (see Section [2.3.21 ). At 
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Table 1 Bulk motions of the local sample. 


UCAC4 Proper Motions PPMXL Proper Motions 


T e g 

(K) 

[Fe/H] 

(dex) 

Number 
of stars 

(Vr) 

(km s -1 ) 

(V*)- 220 
(km s" 1 ) 

(Vz) 

(km s -1 ) 

(Vr) 

(km s" 1 ) 

(Vi) - 220 
(km s -1 ) 

(Vz) 
(km s -1 ) 

(6000, 6800) 

> 

-0.4 

92 

—0.3±1.1 

2.2±0.9 

1.6±0.6 

1.0dbl.3 

2.2±0.8 

2.4±0.6 

(5000, 6000) 

> 

-0.4 

2473 

2.2±0.2 

—5.6±0.2 

—0.2±0.1 

2.0±0.2 

—5.7±0.2 

0.4±0.1 

(4200, 5000) 

> 

-0.4 

2581 

4.3±0.2 

—5.6±0.2 

—1.8±0.1 

4.1±0.2 

—6.2±0.2 

—1.4±0.1 

(6000, 6800) 

< 

-0.4 

46 

—6.6±3.0 

—0.6±2.0 

5.1±0.5 

—4.8±3.4 

—1.7±1.9 

4.5±0.6 

(5000, 6000) 

< 

-0.4 

170 

—3.6±1.5 

—11.9±1.4 

1.4±1.0 

—3.4±1.6 

—14.1±1.4 

1.9±0.9 

(4200, 5000) 

< 

-0.4 

81 

1.3±2.2 

—20.0±1.9 

1.8±1.0 

1.2±2.2 

—19.3±2.1 

1.5±1.2 


this distance limit, a 1 mas yr" 1 error in the proper motion corresponds to a ~ 0.7 km s -1 error in the 
tangential velocity perpendicular to the line of sight. 

The local sample is divided into a metal-rich group ([Fe/H] > —0.4 dex) and a metal-poor group 
([Fe/H] < —0.4 dex). Each of the two groups is further divided into three based on the spectral-types, 
the F-type (6,000 < T eS < 6,800 K), G-type (5,000 < T eS < 6,000 K ) and K-type (4,200 < T eS < 
5,000 K). (Vr), (I4j, ) and (Vz) are estimated by the median velocities of stars of each population, and 
the results can be found in Table 1 and Fig.fl6l 

The differences between UCAC4 and PPMXL results for each population are smaller than 
~ 1 km s -1 , or within the random errors. However, the differences between different populations, ei¬ 
ther in (Vr), ( V,p ) or (Vz), are well beyond the errors, at levels of high significance. The variations of 
azimuthal bulk motions are the most obvious, with (V$) — 220 km s -1 ranging from ~ 2 km s -1 for 
the metal-rich F-type stars to ~ —20 km s -1 for the metal-poor K-type stars. This is what expected, 
considering the larger velocity dispersions and asymmetric drifts of older stars. From F-type to K-type, 
(Vr) increases steadly from ~ 0 to ~ 4 km s -1 for the metal-rich populations, or from ~ —6 km s -1 
to ~ 2 km s -1 for the metal-poor populations. For each spectral-type, the metal-poor stars always have 
lower values of (Vr) than the metal-rich ones. Similarly, (Vz) decreases towards later spectral-types 
amongst the metal-rich populations, from ~ 2 km s -1 for the F-type stars to ~ —2 km s -1 for the 
K-type stars. For the metal-poor populations, (Vz) of the G-type stars are lower than that of the F-type 
by ~ 3 km s -1 , but almost equal to that of the K-type. For each spectral-type, metal-poor stars always 
have larger (Vz) than metal-rich ones. 

5 VELOCITY STRUCTURES IN THE GALACTIC ANTI-CENTER 

Stars along the path from the Sun towards the Galactic anti-center have Galactocentric radial veloci¬ 
ties Vr aligned with their line-of-sight radial velocities V r . With an analysis of the line-of-sight radial 
velocities of ~ 700 red clump stars in the direction of Galactic anti-center, Liu et al. (2012) report a 
“velocity bifurcation” at Galactocentric radii 10-11 kpc and attribute it to a resonance feature of the 
central bar. However, the majority stars of their sample fall just around these radii, it is thus possible 
that the bifurcation may not be a localized feature but rather a global one - it was only detected at these 
radii simply because of the limited distance range probed by the sample. 

To clarify the situation, we make use the anti-center sample defined in Section 12.3.31 consisting 
1,716 GK-, 1,402 L-, 359 A-type dwarfs and 958 RC giants. Their distance distributions can be found 
in Panels (a-d) of Lig. [17] With increasing luminosities, the four groups of stars are able to trace a 
contiguous distance range from the solar vicinity of R ~ 8 kpc out to R ~ 12 kpc. 

Panels (e-i) show the stellar number density distributions of the anti-center sample in the ( R , V r ) 
plane. The distributions exhibit a rich variety of velocity structures at different distance regimes. In the 
first regime, 8.0-8.6 kpc, there are four concentrations of GK-type dwarfs at V r ~ —32, —15, 0 and 
18 km s -1 , respectively. They are apparently related to the known moving groups: the —32 km s -1 
concentration is linked with Wolf 630, the —15 km s “ 1 one with Sirius, the 0 km s -1 one with Coma 
Berenices, and the concentration at 18 km s -1 comes from Pleiades (see e.g. Dehnen 1998; Antoja et 
al. 2012). Moving out to the 8.6-9.0 kpc regime, however, the L-type dwarfs show only two peaks at V r 
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~ —30 and —5 km s _1 , respectively. In the third regime, 9.0-9.5 kpc, there is a single concentration of 
A-type dwarfs at V r ~ —10 kms -1 . Then the concentration moves to V r ~ 0 km s -1 between 9.5-10 kpc 
traced by RC stars. In the regime of 10-11 kpc, the “velocity bifurcation” reported by Liu et al. (2012) 
is clearly visible, appearing as two stellar concentrations at V r ~ 20 and —10 km s , respectively. The 
peaking velocities are slightly lower than theirs, by ~ 6 km s -1 , which possibly originates from the 
measurement uncertainties. Further out to 11-12 kpc, the RC giants show a new triple-peaked structure, 
with peaks located at V r ~ 5 km s -1 , 20 km s -1 and 35 km s -1 , respectively. 

Antoja et al. (2012) have investigated the moving groups in different regions of the disk, showing 
that they drift in the velocity space at different radii. Similarly, the anti-center sample reveals a con¬ 
tiguous spatial variations of velocity structures, from the solar vicinity out to R = 12 kpc. Specially, 
the anti-center sample analyzed here shows that while the “velocity bifurcation” reported by Liu et 
al. (2012) is indeed a localized feature, it is also not the sole feature in the direction of Galactic anti- 
center. It is thus not an obvious conclusion that the “velocity bifurcation” is a resonance feature of the 
central bar, and the formation mechanisms of these features need further investigations. 

6 SUMMARY 

In this paper we have carried out a detailed investigation of bulk motions of the nearby Galactic stellar 
disk using three samples selected from the LSS-GAC DR2: a global sample of ~ 0.57 million FGK 
dwarfs out to a distance of ~ 2 kpc, a local sample of ~ 5,400 FGK dwarfs with distances < 150 pc, 
and an anti-center sample containing ~ 4,400 AFGK dwarfs and RC giants extending out to ~ 4 kpc in 
the Galactic anti-center direction. 

With the global sample, we have provided a three-dimensional map of bulk motions of the Galactic 
disk from the solar vicinity out to ~ 2 kpc with a spatial resolution of ~ 250 pc. Typical values of 
radial and vertical bulk motions range from —15 to 15 km s -1 , while the lag behind the circular speed 
dominates the azimuthal bulk motion by up to ~ 15 km s . The map reveals spatially coherent, kpc- 
scale stellar flows, with typical velocities of a few tens km s -1 . We have analyzed the bending and 
breathing modes and higher-order perturbations. The bending and breathing modes are clearly visible, 
and vary smoothly across the disk plane. The data also reveal breaks and ripples in the profiles of vertical 
bulk motion versus height, indicative of higher-order perturbations. 

With the local sample, we have studied the bulk motions of different stellar populations. Differences 
of several km s -1 are found in the radial and vertical bulk motions, while the azimuthal bulk motions 
exhibit larger differences. 

The anti-center sample reveals abundant velocity structures at distances ranging from the solar 
vicinity out to a Galactocentric radius of 12 kpc, with the nearer ones apparently related to the known 
moving groups. The “velocity bifurcation” reported by Liu et al. is confirmed to be a localized feature 
at 10-11 kpc, beyond which the data reveal a new triple-peaked structure. 
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Fig. 5 Same as Fig.[4]but for the R-Z plane. Numbers in the brackets show the ranges of <1> 
in units of degrees. 
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Fig.6 Same as Fig. |4] but for (V$)—220 km s 1 . 
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Fig. 7 Same as Fig. [5] but for (V$)— 220 km s 1 . 
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Fig.8 Same as Fig. [4] but for (Vz). 
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Fig.9 Same as Fig.[5]but for (VQ). 
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Fig. 10 Distributions of the magnitude and direction of the total velocity \J ( Vr } 2 + ( Vz ) 2 
in the R-Z plane. The arrows indicate the flow directions and their lengths are proportional 
to the total velocities. Red arrows are calculated with the UCAC4 proper motions, whereas 
those in blue are results from the PPMXL. Numbers in the brackets show the ranges of $ in 
units of degrees. The black dot denotes the location of the Sun. 
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Fig. 11 Same as Fig.flOlbut for the total velocity a/(Fr) 2 + ((V$) — V c ) 2 in the X—Y plane. 
Numbers in the brackets show the ranges of Z in units of kpc. 
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Fig. 12 Profiles of ( Vz) as a function of Z. In each panel, the first brackets give the range 
of R in units of kpc, and the second that of $ in units of degrees. Curves with red error bars 
have used the UCAC4 proper motions, while those with blue ones have used the PPMXL 
proper motions. The vertical error bars denote the 1-er random errors of (Vz), whereas the 
horizontal ones are set to the binsize A Z. The UCAC4/PPMXL results have been shifted by 
+/—5 km s -1 , respectively, for clarity. 
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Fig. 13 Breathing mode parameter A in units of km s _1 kpc 1 (colorscale) derived with 
the UCAC4 (lower panel) or the PPMXL (upper panel) proper motions. Each sub-panel plots 
(Vz) as a function of Z (points) along with a linear fit to the data (straight line). The horizontal 
axis of each sub-panel ranges from —2 to 2 kpc and the vertical axis from —30 km s -1 to 
30 km s -1 . Sub-panels enclosed in dotted black boxes show bins where the differences of A 
yielded by the UCAC4 and PPMXL proper motions differ by less than 5 km s -1 kpc -1 . 
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Fig. 14 Same as Fig. [T3]but for bending mode parameter B in units of km s 1 . Sub-panels 
enclosed in dotted black boxes show bins where the UCAC4 and PPMXL values of B from 
UCAC4 and PPMXL differ by less than 5 km s" 1 . 
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Fig. 15 Panel (a): (Vz) as a function of Z for 9.25 < R < 9.5 kpc and —2 < $ < 0°, same as 
for Panel (w) of Fig.|T2] Panel (b): Distribution of g — r colors (blue curves) and r magnitudes 
(red curves) as a function of Z, where the central lines denote the median values of g — r or 
r and 25% stars have values above/below the upper/lower lines. Panels (c, d): Locus of the 
median values of g — r and r in the (g — r, r ) Hess diagram. The colorscale denotes the 
median differences between the UCAC4 and PPMXL proper motions, same as the colorscale 
in Fig. |H 
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Fig. 16 (Vr), (V$) and (Vz) for different populations of the local sample. The red symbols 
are calculated with the UCAC4 proper motions while those in blue with the PPMXL values. 
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Fig. 17 Panels (a-d): Distributions of Galactocentric distances of GK-, F- and A-type dwarfs 
and RC giants in the anti-center sample. The red lines in the four panels mark the distance 
ranges of the boxed regions in Panels (e-h), respectively. Panels (e-h): Stellar number density 
distributions of the four populations in the ( R, V r ) plane. The colorscale is proportional to the 
number of stars in a 500 pc x 10 km s -1 running box relative to the maximum. The white 
lines are the corresponding 40%, 60% and 80% contours from outside to inside. Panel (i): 
Mosaic of the boxed regions in Panels (e-h). The yellow dashed lines and boxes mark the 
boundaries of distance regimes mentioned in Section^ 




































